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We introduce the concept that there are two generic classes of Mott insu- 
lators in nature. They are distinguished by their responses to weak doping. 
Doped charges form cluster (i.e. distribute inhomogeneously) in type I Mott 
insulators while distribute homogeneously in type II Mott insulators. We 
present our opinion on the role inhomogeneity plays in the cuprates. 

1. INTRODUCTION 

Understanding of Mott insulators is one of the central problems of con- 
densed matter physics. In addition to electronic Mott insulators recently 
an interesting Bose Mott insulatorhas been produced by trapping bosonic 
neutral atoms in an optical lattice.£l The recent upsurge in interest in Mott 
insulators is largely due to their dramatic response to doping. For example by 
doping the transition metal oxides such as the cuprates and the manganites, 
one obtains superconductors with the highest known transition temperatures 
and the world's most magnetoresistive materials. 

Mott states, in addition to being insulating, can be characterized by 
the presence or absence of a spontaneously broken symmetry (e.g. spin 
antiferromagnetism) , by the nature of the low energy excitation spectrum 
(e.g. gapped or gapless),! and most recently, bv the presence or absence 
of topological order and charge fractionalization.l3 To this list, we add a 
"type" index that classifies Mott insulators into two types depending on 
their response to doping. 

In a type I Mott insulator an increasing chemical potential induces a 
first order phase transition from an undoped state to a charge-rich state 
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so that the charge density changes discontinuously. In other words, there 
is a range of "forbidden charge density," and two-phase coexistence. In a 
type II Mott insulator, charges go in continuously above a critical chemical 
potential. Depending on the ratio between the delocalization energy and 
the interaction energy these doped charges can form a Wigner crystal or a 
homogeneous liquid phase.Q An important implication of this work is the 
existence of a generic class of Mott insulators which become inhomogeneous 
upon doping. 

For the transition metal oxides discussed above, there exists cpjjsider- 
able evidence that light doping induces spatial inhomogeniety.B&0'Ba Con- 
versely, for Mott insulator such as Sri-xLaxTiOs,iB conventional Fermi liq- 
uid behavior is observed for doping as low as x = 5%, indicative of homo- 
geneity. Thus, the former are type I, the latter is type II. 

Generally speaking, insulating states with charge gap including both 
bancO and Mott insulators occur in crystalline systems at isolated ratio- 
nal "occupation numbers," u = i'*, where i' is the number of particles per 
chemical unit cell. By "doping" we mean a process which causes the occu- 
pation number to shift away from i'* . When lattice translation symmetry 
is not spontaneously broken, v* is typically an integer for bosons, and an 
even integer for fermions with spin. The fermionic state thus may be adia- 
batically connected to a weakly interacting band-insulator, however a Bose 
Mott insulator is always a strong correlation effect. Charge gapped insulat- 
ing states can also occur when z/* is a rational fraction (for fermions this 
includes odd-integer). Usually when that happens, translational symmetry 
is spontaneously broken so that the unit cell of the reduced translation group 
has integral u for bosonal3ll2l and even integral u for fermions. For instance, 
electronic Mott insulators with v* = 1 often exhibit antiferromagnetic N'eel 
long-range order, which doubles the unit cell leading to an effective fg// = 2. 
Nevertheless there exists model bosonic systems for which the Mott state can 
be shown to have no broken symmetries for v* = l/2i3ll30 . (Currently, no 
laboratory system has been found which unambiguously exhibits this exotic 
behavior.) 

The rest of this paper addresses the response of Mott insulators to light 
doping (i.e. i' — > v* — e). Our central observation is that regardless of 
their classification (i.e. symmetry, gap, topology) Mott insulators can be 
divided into two groups differentiated by whether they remain homogeneous 
after doping.ll3 Furthermore we shall show that these two types of insu- 
lating state is analogous to the two types of superconducting state under 
magnetic field. In particular if the Mott insulators are two dimensional and 
the constituent particles are bosons there exists a mathematical mapping, 
the so-called "duality transformation", that relates their zero-temperature 
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response to doping to the finite-temperature response of a 3D supercon- 
ductor to a magnetic field. (The duality transformation has been used 
to deduce two types of doping behavior bv Balents et al in the context of 
doping a spin liquid called "nodal liquid"ll3)The following table summarize 
correspondence between the two. 



T = properties of 2D Bose Mott insulators 


T > properties of 3D Superconductors 


Doping 


Applying magnetic field 


Chemical potential fj, 


Applied magnetic field H 


Induced particle density p 


Magnetic induction B 


World line of doped particles 


Flux tubes 


Quantum delocalization of doped particles 


Thermal meandering of flux tubes 


Type I Mott insulator 


Type I superconductor 


Mott gap 


He 


Effective attraction between doped particles 


Positive N-S interface energy 


Type II Mott insulator 


Type II superconductor 


Effective repulsion between doped particles 


Negative N-S interface energy 


Mott gap 


Hci 


Wigner crystal of doped particles 


Abrikosov flux lattice 


Superfluid state 


Entangled vortex fluid 


Critical /i at which Wigner crystal melts 


Hc2 



2. DOPING A TRANSLATIONALLY INVARIANT BOSE 

MOTT INSULATOR 

In order to avoid the issue of spontaneous symmetry breaking in elec- 
tronic Mott insulators and focus on the absolute essentials, we consider 
the simplest kind of Mott insulators - the ones formed by spin zero point 
bosons on a lattice. Due to Ref.0 this consideration is no longer an academic 
exercise. Consider the following Hamiltonian 

<ij> i 

+\Y.^ij^to'ia^aj -/i^a+Oi, (1) 

where i,j label the lattice sites on a D-dimensional hypercubic lattice, and 
af creates a boson at site i. The first term of Eq. (|l|) describes the quantum 
mechanical "hopping" of bosons from a site i to its nearest neighbors j, 
the second and the third terms describe the pair-wise interaction between 
bosons. The U term is a contact interaction and the Vij terms describe 
interaction between bosons separated by |rj — rj\. 
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Let us focus on the limit where U/t oo. In this "hardcore" hmit 
doubly occupied site costs energy U, hence are excluded. In this case for 
large positive /i there is an unique ground state 

|Mott >= [|a+|0 >, (2) 
j 

in which each site is occupied by one and only one boson and hence u = 1. 
Since this state is separated from all other states with the same particle 
number by an energy gap of order [/, clearly we have an insulator. 

The behavior of the system upon decreasing (i.e. doping) depends on 
the values of Vij. For Vij = the doped holes (i.e the empty sites) interact 
only through the hardcore exclusion. For appropriate ^ where the ground 
state has a small 1 — z^, the hopping term is clearly minimized by delocalizing 
the holes around, since there is plenty of room available for them. In fact it 
is known that such system is an uniform super fluid with superfluid stiffness 
proportional to t{l — v)lnln[l/{l — v)\xJ\ Conversely, if Vij is attractive 
(negative) with range |rj — rj| < i? it is clear that the holes will cluster when 
the strength of this interaction is strong compared with t - the holes will 
phase separate. Thus depending the value of Vij Eq. (0) can describe either 
a type I or type II Mott insulator when v = \. 

Consider, for example, the case in which Vij = —V for {ij) nearest- 
neighbor sites, and Vij = otherwise. Technically in this limit Eq. (|l|) is 
equivalent to the S=l/2 ferromagnetic XXZ model in a z-direction magnetic 
field 

H = -J.y E ^SfS^ + SfS]) - E S^SJ - E (3) 
<ij> <ij> i 

The mapping between these two models relates Jxy to t, to V, hz to fi + 
Vc/2 (c = the coordination number). The z component of the magnetization 
is related to the boson density according to = {v — \/2)N. (Here N is the 
total number of lattice sites.) For < Jxy {i-e.V < t), this model has XY 
order in the absence of ■ In this range of parameters varying causes the 
magnetization to vary continuously. Thus for each fixed the ground 
state is uniform and ferromagnetic. In terms of bosons this means that the 
doped system is a uniform superfluid for all z/. Conversely, for > Jxy (i.e. 
V > t) the model is effectively an Ising ferromagnet with fully polarized 
ground state. In this case exhibits a discontinuity AM^ = N at = 0- 
Thus for all |M^| ^ 1 the ground state exhibits phase separation into two 
oppositely polarized domains. In terms of the bosons, these two domains 
are Mott insulating (zv = 1) and empty (z/ = 0) respectively. 
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3. DUALITY TRANSFORMATION FOR THE BOSE MOTT 

INSUALTOR IN D=2 

While the above discussions applies to arbitrary space dimensions, there 
is a particularly convenient way of thinking about Bose Mott insulators that 
is specific to D = 2. It turns out that for the class of model given by 
Eq. (D, there is a mathematical mapping, the "duality" transformation, 
that provides us an alternative view of the physics of Eq. in terms of the 
vortices of the boson field. It is this mapping that enables us to establish 
a precise connection between the two types of Mott insulators with type I 
and type II superconductors .0113' El In the following we discuss the physical 
content of the duality transformation without going into its technical details 
of the tranformation.EJ 

A vortex is a topological defect in the Bose field. When a boson is adi- 
abatically transported around a vortex, the boson wavefunction acquires a 
phase factor - the Aharonov-Bohm phase. In the dual picture, the vortices 
are viewed as particles (they turn out to have Bose statistics as well), and 
when they are brought around an original boson they acquire an Aharonov- 
Bohm phase. The fact that a boson and a vortex acquire a phase when 
they go around one another implies that bosons and vortices can not Bose 
condense simultaneously. As is well known, in the Bose superfluid phase the 
vortex density fluctuation must be absent. Conversely, in the dual phase, 
where the vortices form a superfluid, the boson density (and hence the dual 
magnetic flux) must be frozen; thevortex superfluid phase is the Mott insu- 
lating phase of the original bosonsBll It is important to note that the absence 
of boson density fluctuation is a necessary but not sufficient condition for 
vortex condensation. For example a static boson density can still frustrate 
vortex condensation because it acts like a background magnetic field. How- 
ever when the static boson density corresponds to an integral v, the vortices 
see a background magnetic flux corresponding to integral number of flux 
quanta per plaquette. (The vortices live on the dual lattice, i.e., the centers 
of the square plaquettes.) This type of flux is "invisible" to the vortices be- 
cause they can be "gauged away" . The u = 1 Bose Mott insulator, discussed 
in the previous section, corresponds to precisely this situation. 

When the boson density is a fraction (i/ = p/q) it is also possible for 
the vortices to condense. Such a state is most naturally accompanied by 
spontaneous translation symmetry breaking, as discussed above, leading to 
an enlarged unit cell with an effective integer uS^ However, it is possible 
to imagine a more exotic Mott state at i' = p/q in which the translation 
symmetry is unbroken. This could happen if q elementary vortices form a 
bound-state, and these composites then condense. Such a composite conden- 
sation is unfrustrated by an unform static boson density v = p/qB Moreover, 
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since the condensate vortex consists of q elementary vortices, the charge 1/q 
bosonic sohton excitation (viewed by the condensate vortices as a flux quan- 
tum) can become a finite energy excitation. Thus there is fractional charge 
solitons! 

In short, a vortex condensate requires the bosons to Mott insulate. Con- 
sequently we can view a boson Mott insulator as a vortex superconductor. 
Doping changes the average background boson density. To the vortices this 
appears as a change in the background magnetic field. In this way the doping 
properties of the boson Mott insulator is related to the magnetic properties 
of the vortex superconductor. 



4. THE ANALOGY WITH SUPERCONDUCTORS IN A 

MAGNETIC FIELD 

In the above discussion a zero temperature boson Mott insulator is 
mapped onto a zero-temperature vortex superconductor (with quantum fluc- 
tuating two-dimensional electromagnetic fields). EHI The final step is to realize 
that the quantum partition function of a (particle-hole symmetric) two di- 
mensional superconductor with fluctuating gauge field is equivalent to the 
classical (i.e. thermal) partition function of a three-dimensional supercon- 
ductor with thermally fluctuating magnetic field. l3 This final correspon- 
dence between the Mott insulator and the classical fiuctuating 3D supercon- 
ductor is summarized in table I. 

A lot is known about thermally fluctuating 3D superconductors. Mean- 
field theory predicted that there are two types. For a type I superconductor 
the magnetic induction B jumps discontinuously from B = to B = H at 
H = He- As the result there is a range of B (i.e. < B < He) in which 
phase separation occurs. For example one way of enforcing a fixed average 
magnetic induction is to place a fiat slab of type I superconductor under 
magnetic field H < He- It is known that when that is done "intermediate 
state" where superconducting region and normal region alternate occurs. 



It is interesting that among many possible in 



lomogeneous structures the 



lammar structure (or stripe) has bee^ observed.E Translate this using table 
I we conclude that a type I Bose Mott insulator undergoes a first order 
insulator — > superfluid transition as a function of chemical potential. The 
density of doped bosons jumps discontinuously at the transition. When the 
doping density is fixed at a value smaller than the critical density at the first 
order transition the system phase separates. 

Mean-field theory predicts that in a type II superconductor the magnetic 
induction B increases from zero continuously at the lower critical field Hd- 
For H > Hci the magnetic induction appear in the form of fiux tubes each 
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enclosing a single quantum of magnetic flux. These flux tubes form a regular 
lattice in the absence of disorder. Unlike the type I superconductors thermal 
meandering of the flux tubes can affect the physics of type II superconductors 
dramatically near Hd and Hc2- For example near Hd the distance between 
neighboring flux tubes is much greater than the range of their interaction 
(the London penetration depth A). As the result thermal meandering of the 
flux tubes can melt the flux lattice into flux liquid .H At larger magnetic 
field the density of flux tube becomes higher so that their interaction can 
stabilize the flux lattice. This flux lattice persists until H — s- Hc2 where 
the thermal meandering melt the flux lattice again. Translating the above 
using table I implies that the extra carriers enter a type II Mott insulator 
continuously at a critical chemical potential. Once they enter they can either 
delocalize (hence Bose condense) or form a Wigner crystal depending on the 
relative importance of the delocalization energy and the interaction energy. 
At at very low carrier density the delocalization always win and we expect 
the boson hopping to render the system a superfluid. 

What determines a superconductor to be type I or type II is the ratio 
between the London penetration depth and the core size of the vortices, 
or more physically in terms of the sign of the interface energy between the 
normal and superconducting regions. Type I superconductors have a positive 
interface energy while type II superconductors have a negative one. As the 
result the flux tubes effectively attract each other in type I superconductors 
while repel each other in type II superconductors. As we have seen this is 
exactly how we turn a type I Mott insulator into a type II one in Eq. ([l|). 

5. DOPING A BOSE MOTT INSULATOR WITH AN ORDER 

PARAMETER 

As we discussed at the beginning of this paper a Mott insulating state 
can be accompanied by translation symmetry breaking. For example let 
us consider Eq. (jfj) with [//t — s- oo and Vij = +V for nearest neighbor 
< ij > and otherwise. For sufficiently strong V and fj, = the ground 
state breaks translation symmetry and bosons form a checkerboard lattice 
and Mott insulate. In this two-fold degenerate ground state the unit cell is 
doubled. Is this Mott insulator type I or type II? 

With the above specific choice of U and Vij Eq. (|l|) is equivalent to 
Eq. dll) with Jxy = —t, Jz = V and h = — Vc/2. For this choice of 
parameters it is known that as a function of hz Eq. (^) exhibits a "spin 
fiop" transition from the antiferromagnetic Ising {Sz) ordered phase into the 
ferromagnetic XY ordered phase. Translate this into the boson language it 
implies that the boson Mott insulator is type I. 
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Before closing this section we argue that the existence of an order param- 
eter (tp) in the insulating state (for the above example ip is the two-sublattice 
density wave order parameter) has an effect in determining the type of the 
Mott insulator. If we assume that doping changes the value of ip, then at 
the interface between doped and undoped region a spatially varying ip is 
necessary. Through the spatial gradient energy 



a positive contribution to the surface energy is resulted. Consequently the 
presence of an order parameter in the insulating state drives the system 
toward a type I Mott insulator. 

6. What ROLE DOES INHOMOGENEITY PLAY IN THE 
CUPRATES; PAIRING IN THE HOLE-RICH ISLANDS IN 

Bi2Sr2CaCu20s+x 

Finally we come to address the question what role does inhomogeneity 
play in the curpates? One of the most direct evidence of electronic inho- 
mogeneity in the cuprates comes from the STM image of the surface of 
Bi2Sr2CaCu20s+x, where nanoscale spatial variation of the energy gap in 
the tunnelling density of states is observed.&i (Of course we have to assume 
that the STM results is not a surface artifact.) As emphasized by Lang 
et alfi the characteristics of the tunnelling spectra divide into two distinct 
types: "a region" where sharp coherence peaks exist (A ranges from 25 -50 
meV), and "/?" region where there is no sharp coherence peak ( A ranges 
from 50 to 75 meV). We take these evidences as suggesting that in the 
a region superconducting pairing has clearly won over all other competing 
orders, while in the f3 region it has not. We believe that inhomogeneity 
allows the superconducting a region to exist in the underdoped regime. As 
the result superconductivity prevails over a wider range of doping concen- 
tration. However, we do not think inhomogeneity plays an important role in 
the pairing of a regions. The above line of thinking suggests that in order to 
understand the pairing mechanism it is best to concentrate on the a regions. 

Since the superconducting gap can vary over the length scale of a few 
nanometers, the coherence length must not be significantly longer than that. 
At such length scale (comparable with the averaged inter-hole distance) the 
Coulomb interaction is poorly screened. How can pairing tolerate such strong 
unscreened Coulomb interaction? 

In the rest of this section we present an explicit example of a paired 
state that can survive strong repulsion. The state we shall concentrate on is 
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t insulator 



a) with projection b) no projection 




Figure 1 

Fig. 1. Results for x — 12%. Ai? as a function of A from Gutzwiller projection 
in a 10 X 10 lattice. Open circles: pure t-J model; crosses: t-J model with nearest 
neighbor repulsion y„„ = 3 J; squares: t-J model with Coulomb interaction Vc = SJ. 

of the form 

\i;{A) >= PimMdBCsi^) > . (5) 

In Eq. (H) \^pdBCs{^) > is a d-wave BCS state characterized by the gap 
parameter A. The operator P is a Jestrow correlator that suppresses con- 
figurations with large density fluctuations. A particular form of P, that is 
often used in the high Tc context, is the Gutzwiller projection operator. In 
our opinion although we do not know precisely what Hamiltonian gives rise 
to a ground state like Eq. it is rather safe to assume that the ground 
state in the a region is of that form. 

In the following we demonstrate that a state like Eq. (^) has the property 
that it is more sensitive to magnetic interactions than to charge interactions. 
To do that we consider the following Hamiltonian 

H = -t^{c^j^Cia+h.c.) + J^{Si-Sj-^inj) 

+ Vc^—{ni-n){nj-n). (6) 

i>j 

The purpose of the following discussions is not to prove Eq. (^) is the ground 
state of Eq. @. Rather we hope to demonstrate one point - the state given 
in Eq. is more sensitive to the magnetic interaction J than to the charge 
interaction Vc- The result we shall quote below is obtained by Wang and 
collaborators in Ref.[26]. 

By a straightforward Monte-Carlo minimization Wang et al conclude 
that for doping x = 0.12 it is energetically favorable to develop a non-zero 
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A for Vc as big as 9J. To appreciate the effects of the Jestrow correlator P 
(which they take as P = YI^{1 — nj-fnj|) in Eq. (|5[) they compared the results 
with and without P. Interestingly without P a nearest-neighbor repulsion 
Vnn destabilizes pairing when Vnn is larger than ~ 0.5J! Thus the Jestrow 
correlator P enhances the stability of pairing in Eq. (^) by as much as a 
factor of 18! 

The following is a brief summary of Wang et aZ's results. Given Eq. (|5|) 
Wang et al minimize E{A) = {^\H\^) / {^\^) by varying A. The results 
presented below are obtained for x = 0.12 in 10 x 10 lattice by variational 
Monte-Carlo. 

In Fig. 1(a) AE = E{A) — E{0) versus A is plotted. The wavefunction 
used is Eq. (^) with P being the Gutzwiller projector. The open circles are 
for the pure t-J model (i.e. without charge-charge interaction), the crosses 
are for t-J model with a nearest neighbor repulsion Vnn = 3 J, and the open 
squares are for t-J -|- Coulomb model (Eq. (^)) with Vc = 3 J. For each of 
the three cases a nonzero A develops. 

In Fig. 1(b) AE = E{A) — E{0) versus A is plotted for the repulsive 
nearest-neighbor (Vnn) model. An important difference is that the Jestrow 
correlator P in Eq. (|5|) is removed. This time the optimal A vanishes for 

Vnn > 0.5J. 

These results clearly indicate that a strongly correlated paired state 
such as Eq. (^) can indeed sustain strong repulsion between the electrons. 
In addition we take this as suggesting a state like Eq. ^ is more sensitive 
to magnetic interactions than the charge interactions. 

In the literature it is often stated that pairing correspond to real space 
binding of holes. In the presence of strong Coulomb interaction real space 
bound hole pairs are extremely energetically unfavorable, and indeed in an 
ansatz such as Eq. (^) no hole binding is present. 

7. CONCLUSIONS 

To summarize, in this paper we introduce the concept that there are 
two types of Mott insulator in nature. A type I Mott insulator becomes 
inhomogeneous after doping, while a type II Mott insulator remains homo- 
geneous. We present a specific lattice boson models which Mott insulates 
without symmetry breaking at v = 1 and, depending on the sign of a mi- 
croscopic interaction, exhibits these two types of ground states after doping. 
We argue that the presence of a symmetry breaking order parameter in the 
insulating state has the effect of driving the system toward type I. In ad- 
dition, we argue that these two types of Mott insulating states are dual to 
the type I and type II superconducting states. Here we conjecture that the 
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notorious x = 0.19 "quantum critical point" in the cuprates is related to the 
inhomogeneous to homogeneous transition when the magnetic induction (B) 
in a type I superconductor is varied across He- 

We also express our opinion on the role played by the inhomogeneity in 
the cuprates. To reiterate, we do not think inhomogeneity plays a central 
role in pairing. However it does allow the superconducting region (where 
pairing dominates over other competing orders) to protrude into the under- 
doped regime. As the result superconductivity prevails over a wider range 
of doping concentration. Finally we believe a correlated pairing state given 
by Eq. (|5|) can describe the hole-rich superconducting regions (the so-called 
a regions) imaged by the STM. Using the result of Ref.[26] we demonstrate 
that this state is much more sensitive to spin-spin interaction than charge- 
charge interaction. This provides an mechanism by which pairing can survive 
poorly screened Coulomb interaction. 

ACKNOWLEDGMENTS 

Part of this paper (two types of Mott insulator) is the result of a collab- 
oration with S.A. Kivelson. I thank Seamus Davis and members of his group 
for valuable discussions. DHL is supported by NSF grant DMR 99-71503. 



REFERENCES 

1. M. Greiner et al, Nature, 415, 39 (2002). 

2. A Mott insulator necessarily has a charge gap, however it is possible for it to 
have a gapless spin excitation spectrum. 

3. T. Senthil and M.P.A. Fisher, Phys. Rev. B 63, 134521 (2001); C. Nayak a nd 
K.Shtengel, Phys. Rev. B 64, 064422 ( 2001); X-G Wen, |cond-mat/010707ll ; L. 



Balents, M. P. A. Fisher, S. M. Girvin, |cond-mat/01 10005 . 

4. Presumably, in the absence of long-range forces, at low enough doping density, 
the charges in a type II Mott insulator always form a uniform fluid, while in 
the presence of long-range Coulomb forces, and at T = 0, they always form a 
Wigner crystal. 

5. For a review of stripes in the cuprates see V. J. Emery, S. A. Kivelson, J. M. 
Tranquada Proc. Natl. Acad. Sci. USA 96, 8814 (1999). 

6. C. Howald, P. Fournier, and A. Kapitulnik, Phys. Rev. B 64, 100504 (2001); S. 
H. Pan et al. Nature, 413, 282 (2001). 

7. K.M. Lang et al, Nature 415, 412 (2002). 

8. For a review of the possible inhomogeneity in the magnites see, e.g., E. Dagotto, 
T. Hotta, A. Moreo, Physics Report 344, 1 (2001). 

9. See, e.g., J.M. Tranquada, J. Phys. Chem. Solids, 59, 2150 (1988). 

10. Y. Tokura et al, Phys. Rev. Lett. 70, 2126 (1993). 

11. A band insulator insulates because of Pauli exclusion principle. In the non- 
interacting prototype of this type of insulator all available single-particle orbitals 



D.-H. Lee 



are occupied. This type of insulator usually possesses a even number of electrons 
per crystalline unit cell. 

12. D-H Lee and R. Shankar, Phys. Rev. Lett., 65, 1490 (1990). 

13. M. Oshikawa, Phys. Rev. Lett. 84, 1535 (2000). 

14. R. Moessner and S.L. Sondhi, Phys. Rev. Lett. 86, 1881 (2001). 

15. D.S. Rokhsar and S.A. Kivelson, Phys. Rev. Lett. 61, 2376-9 (1988). 

16. In the presence of long-range Coulomb forces macroscopic charge inhomo- 
geneities is impossible as — z^* — > 0. However, so long as the long-range forces 
are not too strong, there could exist a range of low doping in which the doped 
charges tend to cluster. We still refer to this kind of system as type I Mott 
insulators as well. For instance, if a low density of doped charges form puddles 
or stripes, we consider this a form of micro-phase-separation, while a Wigner 
crystal of doped charges is deemed homogeneous. This is because in a Wigner 
crystal charges localize in the smallest unit hence no clustering is exhibited. 

17. D.S. Fisher and P.C. Hohenberg, Phys. Rev. B 37, 4936 (1988). 

18. M.P.A. Fisher and D-H Lee, Phys. Rev. B 39, 2756 (1989). 

19. L. Balents, M.P.A. Fisher and C. Nayak, Int. J. Mod. Phys. B 12, 1033 (1998). 

20. The system considered in this Ref.El can be "doped" by varying the period of 
the trap lattice. 

21. We note that the vortex superfluid discussed here differs from the generic boson 
superfluid because the net vorticity is zero. The field theory that describes such 
" particle- hole symmetric" condensate has space and time on equal footing (i.e. 
relativistic) . This fact makes the later analogy with thermal fluctuating 3D 
superconductors completely appropriate. 

22. The field strength F,,^ (^ = 0,1,2) of a 2-f 1-dimensional EM field has three 
independent components. The correspondence of the these three components 
with the boson three-current (p, j) is as follows: F12 p, Fqi ^ j2, -F02 ji- 

23. In the correspodence between a quantum superconductor at T = and a ther- 
mally fluctutating superconductor at finite temperature the magnetic field (-^12) 
and electric fields {Fqi, -F02) in the 2-1-1 dimensional quantum theory is mapped 
on to Bz and By, —B^ respectively. 

24. See. e.g., M. Tinkham "Introduction to superconductivity", McGraw-Hill, New 
York (1975); P.-G. de Gennes, "Superconductivity of metals and alloys", Ben- 
jamin, New York (1966). 

25. D.R. Nelson, Phys. Rev. Lett. 60, 1973 (1988). 

26. Q-H Wang, J. H. Han and D-H Lee, Phys. Rev. B 65, 054501 (2002). 



